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ABSTRACT 

We present detailed morphological properties of Lya emitters (LAEs) at z ~ 
5.7 in the COSMOS field, based on Hubble Space Telescope Advanced Camera 
for Surveys (ACS) data. The ACS imaging in the F814W filter covered 85 LAEs 
of the 119 LAEs identified in the full two square degree field, and 47 LAEs of 
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them are detected in the ACS images. Nearly half of them are spatially extended 
with a size larger than 0.15 arcsec (~0.88 kpc at z — 5.7) up to 0.4 arcsec (~2.5 
kpc at z — 5.7). The others are nearly unresolved compact objects. Two LAEs 
show double-component structures, indicating interaction or merging of building 
components to form more massive galaxies. By stacking the ACS images of all 
the detected sources, we obtain a Sersic parameter of n ~ 0.7 with a half-light 
radius of 0.13 arcsec (0.76 kpc), suggesting that the majority of ACS detected 
LAEs have not spheroidal-like but disk-like or irregular light profiles. Comparing 
ACS F814W magnitudes (hu) with Subaru/Suprime-Cam magnitudes in the 
NB816, i', and z' bands, we find that the ACS imaging in the F814W band 
mainly probes UV continuum rather than Lya line emission. UV continuum 
sizes tend to be larger for LAEs with larger Lya emission regions as traced by the 
iV.B816 imaging. The non-detection of 38 LAEs in the ACS images is likely due 
to the fact that their surface brightness is even too low both in the UV continuum 
and Lya emission. Estimating Igu for the LAEs with ACS non-detection from 
the z' and iV.B816 magnitudes, we find that 16 of these are probably LAEs with 
a size larger than 0.15 arcsec in UV continuum. All these results suggest that 
our LAE sample contains systematically larger LAEs in UV continuum size than 
those previously studied at z ~ 6. 

Subject headings: cosmology: observations — cosmology: early universe — galax- 
ies: formation — galaxies: evolution — galaxies: morphology 

1. INTRODUCTION 

During the last decade, a large number of young, star-forming galaxies beyond redshift 
of z = 5 have been found based on deep imaging observations of both the Hubble Space 
Telescope (HST) and 8- 10m class telescopes such as the 8.2m Subaru Telescope (see for a 
review, Bouwens & Illingworth 2006; Taniguchi 2008; see also for recent progress, Bouwens et 
al. 2008, 2009; Bradley et al. 2008). In particular, narrow-band imaging surveys have been 
providing us well-defined samples of strong Lya emitters (LAEs) at z ~ 5.7 (e.g., Rhodes 
& Malhotra 2001; Ajiki et al. 2003; Hu et al. 2004; Shimasaku et al. 2006; Murayama et 
al. 2007 and references therein). These surveys are used to investigate the star formation 
activity in such young galaxies, providing typical star formation rates from several to a few 
tens M yr^ 1 . Clustering properties are also one of the important issues studied as structure 
formation in the early universe provides an important observational constraint on hierarchical 
structure formation scenarios (e.g., Springel et al. 2005). Although a possible over-density 
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region has been identified in the Subaru XMM-Newton Deep Survey Field (Ouchi et al. 
2005), there is no other significant evidence for clustering of young galaxies at z ~ 5.7 in 
other deep wide-field surveys (e.g., Shimasaku et al. 2006; Murayama et al. 2007). 

Another interesting question is addressed to morphological properties of LAEs because 
these properties give us insights on how LAEs were assembled and how their intense star 
formation events were triggered. However, no systematic investigation of the detailed mor- 
phology of high-z LAEs has yet been undertaken although some case studies have been 
reported (e.g., Rhoads et al. 2005; Venemans et al. 2005; Pizkal et al. 2007; Overzier et al. 
2008 and references therein). Rhoads et al. (2005) found in the Hubble Ultra Deep Field 
(UDF) a LAE at z = 5.4 (UDF 5225) with a spatial extent of 0.25 arcsec x 1.0 arcsec (1.6 
kpc x 6.3 kpc). This LAE shows a core together with three knots that appear to emanate 
from the core. On the other hand, Bunker et al. (2003) found a very compact LAE at 
z = 5.78 in the Chandra Deep Field South among their i-dropout sample. Its half light 
radius is only 0.08 arcsec (490 pc), barely resolved by the ACS PSF (0.06 arcsec). Stanway 
et al. (2004a) found three LAEs in the UDF during the course of their GLARE project (= 
Gemini Lyman Alpha at Reionization Era) and obtained also small sizes ranging from 0.09 
arcsec (~ 500 pc) to 0.14 arcsec (~ 900 pc). Relative small sizes of six LAEs with respect to 
Lyman break galaxies (LBGs) at similar redshift (z ~ 6) are also found by Dow-Hygelund 
et al. (2007). Since the size and morphological properties provide us with useful insight on 
the understanding of physical processes of star formation, we need systematic studies of the 
detailed morphological properties for a large sample of such high- z galaxies. 

Recently, Murayama et al. (2007; hereafter M07) identified a total of 119 LAEs at z ~ 
5.7 in the Cosmic Evolution Survey (COSMOS) field (Scoville et al. 2007a), providing one 
of the largest samples of LAEs in a large contiguous field. Since F814W imaging taken with 
the Advanced Camera for Surveys (ACS) on-board the HST is available for the COSMOS 
field (Scoville et al. 2007b; Koekemoer et al. 2007), the sizes and morphologies of the LAEs 
in the COSMOS field can be investigated in detail. In particular, since the redshifted Lya 
emission of our LAE sample is probed in the F814W imaging, it is possible to study the Lyct 
morphology as well as the rest-frame ultraviolet continuum shape. In this paper, we present 
our detailed analysis of ACS images of the LAE sample of M07. 

We use a standard cosmology with f2 ma ttcr = 0-3) = 0.7, and H = 70 km s _1 Mpc -1 . 
Throughout this paper, we use magnitudes in the AB system. 
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2. OBSERVATIONAL DATA AND ACS COUNTERPARTS OF LAEs 

In M07, 119 LAE candidates at 5.65 < z < 5.75 were carefully selected from optical 
imaging with both the narrow-band filter, NB816 (A c = 8150 A with a width A A = 120 
A see Ajiki et al. 2003 for details) and broad-band filters from B to z' taken for a 1.95 
deg 2 area of the COSMOS field using Suprime-Cam (Miyazaki et al. 2002) on the Subaru 
Telescope (Kaifu et al. 2000; lye et al. 2004). Details of the Subaru observations and data 
processing are described by Taniguchi et al. (2007) and Capak et al. (2007). Note that 
follow-up spectroscopy has been performed for 24 LAEs in the sample and that all of them 
showed Lya emission at z ~ 5.7 (Capak et al. 2009) verifying the effectiveness of the adopted 
selection method. 

The ACS data in the F814W filter were taken for an area of 1.64 deg 2 of the COSMOS 
field and were processed to 0.05" pixel -1 images with an averaged point spread function 
(PSF) width of 0.097" (Scoville et al. 2007b; Koekemoer et al. 2007). In our analysis, we 
use the official COSMOS ACS image, Version 1.3. 

Since the observed area of the ACS imaging is slightly smaller than that of the Subaru 
one, we find that ACS data are available for 85 LAEs in the LAE sample defined by M07. 
Spatial distribution of these LAEs is shown with red symbols in Fig. 1. Among them, 20 
LAEs are already spectroscopically identified as star-forming galaxies at z ~ 5.7 (Capak et 
al. 2009). 

Our data analysis procedures for ACS data are as follow. First, for each LAE, we created 
a small ACS cut-out (16"xl6") centered on the the LAE position (as derived in the NB816 
image; see M07). Then, we used SExtractor (Bertin & Arnouts 1996) to detect a possible 
ACS counterpart of the LAE. We adopted the following SExtractor parameters: a detection 
threshold of 1.6er and a minimum detection area of 9 pixels. These values were carefully 
selected (1) to allow for the detection both diffuse elongated sources as well as point-like 
sources, (2) to minimize false detections, and (3) to minimize detection failures for sources 
that are apparently shown in the ACS image by eyes. A gaussian-profile filter (5 pixels x 
5 pixels) with a FWHM of 2 pixels (0.1") well matched to the PSF size was applied to the 
detection images for smoothing. The SExtractor parameters for the source detection in Isu 
are listed in Table 1. 

We found 58 sources detected near the LAE positions (r < 1"). By eye inspection we 
rejected 3 sources because they are largely offset (> 0.68") from the LAE center and have 
no distinct counterpart in the Subaru broad-band images of any bands. Nest, in order to 
remove low- z foreground neighbors, we examined the Subaru B, V, and g' images. Since 
the observed wavelength of the Lyman limit at z ~ 5.7 shifts to 6110 A, a true LAE must 



- 6- 



be undetectable in these bands. In this analysis, we adopted the search radius for ACS 
counterparts is 1" and then removed 6 sources. 

Finally, among the 85 LAEs whose ACS data are available, we found 49 ACS sources 
for 47 LAEs (17 objects were spectroscopically confirmed; Capak et al. 2009). Two LAEs 
(No. 60 and No. 110) have double-component ACS sources. Offsets of the ACS positions 
from the NB816 positions are typically less than 0.13" in the ACS images, smaller than the 
pixel scale (0.15" pixeE 1 ) of the NB816 images. The remaining 38 LAEs (3 objects were 
spectroscopically confirmed; Capak et al. 2009) are not detected in the ACS images. The 
numbers of the total sample, both the detected and not-detected LAEs are summarized in 
Table 2. 

In Fig. 2, we show thumbnails of the 85 LAEs in the ACS F814W images. We also 
show thumbnails in the ACS F814W images together with their Subaru i', NB816, and z 1 
images in Fig. 3. Smoothed F814W images with a gaussian kernel with a FWHM of 2 pixels 
are also presented. The detected sources identified as LAE counterparts are indicated by 
red ellipses on the smoothed ACS, NB 816, and ^'images. Green and blue ellipses are ACS 
sources excluded from the sample by eye inspection and by rejection of foreground neighbors, 
respectively. 

As shown in Fig. 1, the ACS detected (red filled circles) and not-detected LAEs (red 
crosses) appear to be almost randomly distributed in the COSMOS field and thus their dis- 
tributions may not be affected by large-scale inhomogeneity of the ACS data quality (e.g., 
edges of the field). The total magnitude (Jgu) and half-light radius (Rhl) were measured 
for each detected source by SExtractor on the original ACS image (i.e., not on the smoothed 
image). We adopted SExtractor's MAG_AUTO as total magnitude. For the 38 LAEs un- 
detected in the F814W image, we estimated 3a upper-limits for the magnitudes within a 
1" diameter aperture. These photometric properties of the ACS data are listed in Table 3. 
Note that the 3a limiting magnitude of the F814W images is 27.3 mag in a 1" diameter 
aperture. All magnitudes are corrected for the Galactic extinction of A F8UW =0.035 (Capak 
et al. 2007). In Table 4, we list the photometric properties of the LAE candidates from M07. 
The 3a limiting magnitudes within a 2" diameter aperture in the NB816, i', and z' images 
are 25.7, 26.1, and 25.3, respectively. 
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3. MORPHOLOGICAL PROPERTIES 

3.1. Half Light Radius 

The ACS counterparts of the LAEs look differently from object to object as shown in 
Fig. 2. Some LAEs show compact, round, and nearly unresolved shapes, while others have 
an elongated, filamentary, or patchy morphology. 

First, we analyze the sizes of our LAE sample in the ACS images. We show the distri- 
bution of half-light radius (Rhl) for the 47 ACS-detected LAEs in Fig. 4. In this histogram, 
the size of the larger component is adopted for the two LAEs with double-component ACS 
sources. Note that the measured half-light radii of stars in the ACS images are typically 
0.11" as indicated in Fig. 4. The number of LAEs in each /?hl bin decreases with increasing 
half-light radius beyond 0.15" while the overall distribution appears to be almost flat at 
the radius of R^l < 0.15". Therefore, we may conclude that LAEs with i?HL < 0.15" are 
almost unresolved compact objects. Three LAEs (Nos. 40, 68, and 78) have a larger size 
as i?HL > 0.3" compared with the PSF size. However, the sizes of all the LAEs detected in 
our ACS imaging are smaller than 0.4" and no widely extended LAE (like the LAE found 
by Rhoads et al. 2005) has been found. 

In Fig. 5, we show the distribution of Igu- This distribution is affected by detection 
incompleteness toward fainter magnitudes because the detection is limited by surface bright- 
ness. For the LAEs with no detection in ACS, we show 3a upper limits within a 1" diameter 
aperture in Fig. 6. If we adopt la and 2a upper limits, the limiting magnitudes are fainter 
by 1.2 and 0.45 magnitude, respectively. 

In order to examine the effect of limiting surface brightness in our ACS imaging, we 
show the relation between i? H L and J 814 magnitude in Fig. 7. In this diagram, we find no 
faint object with a large i? H L- This suggests that we are unable to detect LAEs with i? H L > 
0.5 arcsec if they are fainter than 26 mag in / 8 i 4 . 

In Fig. 8, the frequency distributions of sizes in NB816 imaging, FWH M(NB816), 
are shown for the detected and not-detected LAEs. Not a small part of the sample LAEs 
have larger sizes than the PSF size (0.98 arcsec), especially for the not-detected LAEs. Their 
appearance in the NB816 images shows widely extended morphology (e.g., No. 74 in the ACS 
detected subsample, and No. 1, No. 27, and No. 29 in the ACS not-detected subsample) and 
looks like a so-called Lyct blob (e.g., Fynbo et al. 1999; Steidel et al. 2000; Matsuda et al. 
2004; Saito et al. 2006). In Fig. 9, we also show the diagram between FWH M(NB816) and 
i?HL- This tendency is more conspicuous for the not-detected LAEs. Even if a LAE has a 
larger size in NB816, i.e., it has an extended Lyct source, we would be able to detect it in 
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our ACS F814W imaging if it has a compact UV continuum component less than 0.5 arcsec. 
Therefore, it is suggested that the not-detected LAEs have an extended UV continuum 
component larger than 0.5 arcsec. However, the z' magnitudes of the not-detected LAEs are 
typically fainter relative to the ones detected in ACS (see Table 3). This means that their 
UV magnitudes are in some cases too faint to be detected in our ACS imaging even though 
they are not very much extended (i.e., < 0.5 arcsec); see section 3.4 for more details. 

3.2. Monte Carlo Simulations 

To estimate practical errors in our measurements of -Rhl, we performed Monte Carlo 
simulations. We prepared 1000 artificial sources with the exponential light profile for each set 
of given parameters of the total magnitude, the half light radius, and the ellipticity. We put 
these sources on the observed ACS image after they were convolved with the PSF image. We 
also added photon noises. Then, we measured their photometric properties with the same 
SExtractor parameters for the source detection. Detection completeness of the 50% limit is 
indicated in Fig. 7. 

Based on these simulations, we estimated probability distributions of each parameter for 
each LAE. The median values with the 68% confidence intervals of estimated distributions 
for total magnitude and half light radius are listed in Table 3. In the upper panel of Fig. 
10, we show the relation between estimated magnitudes and measured magnitudes. The 
estimated magnitudes tend to be brighter than the measured magnitudes. The estimated 
errors are typically ±0.3 - ±0.5 and larger than the measured error (typically ±0.1-±0.2). 
Although the estimated half light radius suffers from large uncertainty, the LAEs measured 
with i?HL > 0.3" still have a large estimated value in -Rhl- Moreover, the simulation indicates 
several LAEs are intrinsically large even their measured sizes are nearly the PSF size. While 
most of LAEs are unresolved objects, we may conclude that some of LAEs at z ~ 5.7 have 
larger intrinsic _R H l than the PSF size. 

3.3. What Do We See in ACS F814W Images; Lya Emission or Ultraviolet 

Continuum ? 

As we see in Section 3.1, the 47 LAEs detected in our ACS F814W imaging are not 
largely extended (i.e., /?hl < 0.4 arcsec). What do we see exactly in the ACS images? Here 
we examine whether the detected light is from Lya emission or Ultraviolet (UV) continuum. 

First, we present the transmission curves for the filters used in our analysis; F814W for 



- 9- 



ACS/HST and i', NB816, and z' for Suprime-Cam/Subaru in Fig. 11. The CCD sensitivity 
is taken into account for each filter transmission curve . The bandpass of the NB816 filter is 
almost centered in the wavelength range cover the F814W filter. Therefore, if Lya emission 
is strong enough to be detected in our F814W imaging, we would see Lya morphologies of 
LAEs. However, the UV continuum at wavelengths longer than 1216 A can also be probed 
by the F814W imaging whose transmission curve is similar to the sum of Suprime-Cam i' + 
z' filter transmission. In order to investigate what our F814W imaging probes, we examine 
the correlation between J 814 and NB81Q in Fig. 12. Since the correlation appears to be poor 
(its correlation coefficient is r = 0.40), the F814W imaging does not primarily probe Lya 
emission. Next, in Fig. 13, we show the correlations between J 814 and i' (r = 0.32), and 
^814 and z' (r = 0.45). These comparisons show that J 8 i 4 is more correlated with z'. This 
suggests that the F814W imaging probes UV continuum from massive stars in each LAE, 
because z' is only sensitive to UV continuum as evident from the filter curves (Fig. 11). 

Although the correlation between 7gi4 and z' is better than for NB816 and i', there is a 
systematic offset between I 8 u and z' values. Isu is typically 0.94 mag (average) fainter than 
z' . This offset can be explained by the difference of their filter transmission curves between 
J 814 and z' (see Fig. 11). Since the Lya wavelength, 1216 A, is observed nearly at the band 
center of the F814W filter, the UV continuum at wavelengths longer than 1216 A is covered 
by only half of the F814W bandpass (see Fig. 11). As a result, flux densities of the UV 
continuum are underestimated with our F814W imaging by roughly 0.75 (=2.5 log 2) mag. 
For a more precise treatment, we can estimate the correction factor that converts J 8 i 4 to z' 
for each LAE by assuming a simple model spectrum similar to the one shown in Fig. 11, 
in which f v is assumed to be constant. Given both the UV continuum flux estimated by 
our z' magnitude and the rest-frame EW(Lja), we estimate the correction factor for each 
LAE and then estimate z 1 magnitude based on J 814 value, z'(/ 814 ). For the LAEs with only 
lower limits of EW (Lya), we use the lower limit values for this estimate. In Fig. 14, we 
compare z'(I 8 u) with z', and find that there is a tighter correlation (r = 0.58) between these 
two magnitudes. Therefore, we conclude that we see the UV continuum of LAEs and miss 
almost all the Lya flux in our ACS imaging. 



3.4. Spatial sizes of the LAEs without ACS counterparts 

It is reminded that 38 LAEs are not detected in our ACS imaging. This means that the 
surface brightness of these LAEs is too low to be detected not only in Lya but also in their 
UV continuum. Even for LAEs with bright UV continuum, they could not be detected in 
the ACS image if they are spatially extended and thus their surface brightness falls below 
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the detection limit. Therefore, some of the LAEs undetected in the ACS images may have 
large .Rhl- 

Here, assuming that the LAEs without ACS counterparts also have the same SED as 
that adopted in the previous subsection, we estimate the expected Isu magnitude, / 814 (V), 
based on the Subaru z' and the iV.B816 magnitudes. We find that 16 of the 38 not-detected 
LAEs in ACS would be bright as much as Isu{z') < 26.5. For LAEs with Jgi4=26.5 and 
i?HL=0.15 arcsec, the detection completeness is 50% (see Fig. 7). Therefore, these 16 LAEs 
with Isi4(z') < 26.5 are probably also slightly extended LAEs (Rul > 0.15 arcsec). However, 
there is no LAE without ACS counterparts brighter than 25 mag in Isia(z'). This suggests 
that very extended LAEs with /?hl > 0.4 arcsec are not present in our LAE sample. 

As shown in the previous subsection, z' is well correlated with J 814 and these magnitudes 
represent UV continuum brightness for the ACS-detected LAEs. For the not-detected LAEs, 
it is not uncertain if this relation would be still valid. Their z 1 are systematically fainter 
than those of the ACS detected LAEs while the magnitude ranges of AS816 for both the 
subsamples are nearly the same, implying that contribution of Lyct flux would be considerable 
for that of the not-detected LAEs. Since the sizes in NB81Q of the not-detected LAEs are 
systematically larger than those of the ACS-detected LAEs (Fig. 8), it is suggested that a 
faint compact UV continuum source with a widely extended bright Lya nebula explains for 
the undetection in /gi4 for a part of the not-detected LAEs. However, at least for the 16 
not-detected LAEs with Isu(z') < 26.5, their rest EW(Lya) is comparable with that of ACS- 
detected LAEs, so that they would be detected in ACS if they had a compact (Rul < 0.15 
arcsec) UV continuum source. Therefore, they probably have an extended UV continuum 
source. 



4. DISCUSSION 

4.1. ACS Size and Star Formation Properties 

In the previous section, we have demonstrated that our ACS imaging with F814W does 
not probe Lya emission but UV continuum from massive stars in the LAEs. Since all LAEs 
detected in the ACS imaging are spatially small (i.e., < 0.4 arcsec or <2.5 kpc), their star- 
forming regions are considered to be physically compact. The sizes of high-z galaxies provide 
important information on the growth of the luminous parts of galaxies embedded in dark 
matter haloes (e.g., Dalcanton, Spergel, & Summers 1997; Mo, Mao, & White 1998; Ferguson 
et al. 2004; Bouwens et al. 2006). 

Here we investigate how the size of the UV continuum is related to the global star 
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formation properties in the LAEs. In Fig. 15, we show the diagram between i?HL and 
L(Lya). In Fig. 16, we present the diagram between i? H L and the rest-frame EW(Lya). 
We find no correlation in these two diagrams. These results suggest that the size of UV 
continuum regions is not directly related to the star-forming activity traced by Lja emission 
in our LAEs. However, one would naturally have expected positive correlations in these 
diagrams. Therefore, the little correlation implies the following three possibilities. (1) Most 
star-forming regions have lower surface brightness than our detection limit. (2) Most star- 
forming regions are hidden by dusty clouds. Or, (3) the UV continuum probes star forming 
regions with a large age spread while Lya only probes youngest star forming regions that 
might not spatially overlap with older star forming regions. Radio and millimeter stacking 
analysis for our LAE sample rules out large dust content (Carilli et al. 2007). Since a number 
of observations also suggest that LAEs at high redshift tend to have little dust content (Lai 
et al. 2007; Finkelstein et al. 2008, 2009a, & 2009b; see, however, Chary et al. 2005), the 
second possibility appears to be unlikely. Currently, however, we have no firm answer on 
this issue. In future, we need deeper F814 imaging of our LAEs. Also, sensitive rest-frame, 
mid- and far-infrared imaging will be necessary to solve this problem. 



4.2. Structural Properties of the LAEs Detected in our ACS Imaging 

As we have shown, our ACS imaging with F814W probes UV continuum from massive 
stars in the LAEs and the majority of our LAEs show extended morphology in J§i4. It is 
of great interest to analyze not only their sizes but also their surface brightness profiles. 
Unfortunately, our ACS images are not deep enough to perform such an analysis for each 
LAE individually. Therefore, we alternatively adopt the image stacking method to obtain 
mean properties for our LAEs at z « 5.7. 

First, we rotated the major axis to align in the x-axis for each LAE image. In this 
procedure, we used the position angles measured by SExtractor. We also constructed second 
rotated images in 180 degree from the first rotated (major- axis aligned) images. Next, we 
generated a composite image by co-adding counts of both the first and second rotated images 
of LAEs. In this analysis, we did not use the double-component LAEs, No. 60 and No. 110. 
We also excluded No. 73 and No. 81 because foreground galaxies are very close to these 
LAE sources. Therefore, we used the remaining 43 LAEs for the stacking analysis. 

Finally, we measured structural parameters on the resultant image by modeling 1- 
dimensional surface brightness profile along the major axis. The PSF image required for 
convolving the model light distribution was derived by combining stellar images near the 
LAEs. Assuming a Sersic function, we obtained the PSF-deconvolved effective radius of 
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Rhl = 0.13l£o? arcsec (0.76 kpc) and the Sersic index (n) of n = 0.7lg;|. In Fig. 17, we 
show the composite images, light profiles of the composite image along the major axis, and 
the model light profile. Note that the observed effective radius of stars are typically 0.11 
arcsec, and thus the composite LAE image is only slightly resolved. The derived parameters 
may contain larger uncertainties than the nominal errors estimated from the residuals. How- 
ever, the derived small Sersic index suggests that our LAEs may have on average irregular 
or disk-like morphologies rather than spheroidal structures. 

Ravindranath et al. (2006) found 40% of Lyman break galaxies (LBGs) at z > 2.5 have 
exponential profiles, 30% of LBGs have r 1//4 -like profiles, and 30% of LBGs have multiple 
cores. In contrast, our analysis indicates that spheroidal structures are rare in our LAE 
sample at z ~ 5.7. Note that clumpy- or chain-like structures are often seen in galaxies 
at z=l-5 in the UDF sample (e.g., Elmegreen et al. 2007). Since our stacking analysis 
smears out such fine structures, we cannot rule out that our LAE sample includes LAEs 
with irregular morphologies. 

4.3. Relationships between LAEs and LBGs 

Since LAEs and LBGs are two major populations of star-forming galaxies at high red- 
shift, it has been often discussed what their physical and evolutionary relations are. LAEs 
are selected by the narrow-band imaging technique that probes strong Lja line emission 
(see for reviews, Taniguchi 2005, 2008). This technique does not require that the UV con- 
tinuum emission is strong enough to be detected in optical broad-band imaging. The larger 
EW(Lja) suggest that they tend be younger in context of the elapsed time from the onset 
of star formation activity (e.g., Rhoads & Malhotra 2001; Malhotra & Rhoads 2002; Nagao 
et al. 2007). Therefore, it is expected that LAEs tend to be young less massive star- forming 
galaxies. On the other hand, LBGs are selected as so-called dropout objects in broad-band 
images (e.g., Steidel et al. 1999). This technique requires that the UV continuum is strong 
enough to be detected. However, the selection of dropouts is not affected by the strength of 
the Lyct line in principle (e.g., Steidel et al. 2000; Shapley et al. 2003). Therefore, LBGs 
tend to be more massive and relatively older than LAEs. 

In fact, several studies of LAEs and/or LBGs support the above difference between the 
two populations at z ~ 3 - 5 as follows. For example, LAEs tend to have bluer UV continua 
than LBGs and this property cannot be explained only with a difference in reddening by 
dust (e.g., Gronwall et al. 2007; Ouchi et al. 2008). However, the number fraction of 
LAE/LBG tends to increase with increasing redshift; i.e., from several to 10% at z ~ 3 to 
~ 30% at z ~ 6 - 7 (e.g., Shimasaku et al. 2006; Dow-Hygelund et al. 2007; Ouchi et al. 
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2008; Sumiya et al. 2009). Since, at z ~ 6, the elapsed time from the Big Bang is at more 
1 Gyr. LBGs and LAEs may tend to share nearly the same physical properties although 
their Lya emission luminosities are different on the average. Therefore, it seems important 
to compare observational properties between LBGs and LAEs at z ~ 6 in a more systematic 
way (e.g., Dow-Hygelund et al. 2007; Pentericci et al. 2007). Motivated by this, we discuss 
the size-mass relation both for LBGs and LAEs using our own data presented here together 
with all available data from the literature. 

Recently, Dow-Hygelund et al. (2007) made a spectroscopic study of 22 LBGs around 
z ~ 6 selected from the samples in three different sky areas and they identified six LBGs with 
Lya emission (i.e., LAEs) at z = 5.5 - 6.1. They investigated the size-magnitude relation for 
both LAEs and LBGs compiled from the literature as well as their own data. They defined 
the size as the half-light radius (Rrl) and the magnitude is ACS 2 8 50 magnitude (see Fig. 
13 in their paper). They found that the LAEs tend to be more compact than the LBGs at 
z ~ 6, suggesting that the LAEs are younger than the LBGs. This interpretation may be 
supported by other observational differences between LAEs and LBGs; LAEs tend to be less 
massive (e.g., Overzier et al. 2006; Gawiser et al. 2006; Pentericci et al. 2007). 

In order to increase the sample size for an analysis of the size-magnitude relation com- 
pared to Dow-Hygelund et al. (2007), we have also compiled available data from the literature 
for z ~ 6. Our data compilations is summarized in Table 5. Since our ACS magnitude is 
not Zg5o but 7814) we have to convert our J 8 i4 magnitude to z^q. In this procedure, we 
assume that the flux of the rest-frame continuum at A < 912A is zero while that at A > 
912A is constant. We also calculate the contribution of the Lya flux by using the rest-frame 
EW(Lya). In this way, we derive the following equation for conversion; 

z 850 = / 814 + 2.51og(0.364 + 2.088 x 10~ 3 £;Wo). 

In Fig. 18, we show our own results together with the compiled data for z ~ 6. Our 
new data presented in the right panel add more information on this size- magnitude diagram 
with respect to that presented by Dow-Hygelund et al. (2007). We find that there is little 
difference between the two populations, LAEs and LBGs at z ~ 6 although there many faint 
LBGs taken from the very deep ACS imaging of LBGs made by Bouwens et al. (2006) and 
Bunker et al. (2003) are located at z 8 50 < 29.5. 

5. CONCLUSIONS 

We have derived the detailed morphological properties of Lya emitters (LAEs) at z ~ 5.7 
in the COSMOS field, based on the HST ACS imaging in the F814W filter. Our main results 
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and conclusions are summarized below. 

(1) Among the 119 LAEs at z = 5.7 identified in the HST COSMOS field (M07), 85 
LAEs are imaged with ACS/F814W. Of those, 47 LAEs are detected in our ACS imaging 
while the remaining 38 ones are not detected. 

(2) All LAEs detected in ACS have small spatial sizes (.Rhl ^0.4 arcsec). However, 
nearly half of them show a spatially extended morphology with effective radii larger than 
0.15 arcsec (0.93 kpc), which is larger than the PSF size (-Rhl = 0.11 arcsec). 

(3) Among the 38 ACS not-detected LAEs, 16 LAEs may be spatially extended (-Rhl > 
0.15 arcsec) as estimated from their z' and A .8816 magnitudes. 

(4) Comparing the ACS data with our Subaru A5816, i', and z' data, we find that the 
ACS/F814 imaging probes not Lyct line emission but UV continuum arising from wavelengths 
longer than 1216 A. 

(5) We find a tendency that LAEs with a larger UV continuum source have a larger Lya 
size as probed by our Subaru NB816 imaging. Since the LAEs with ACS non-detections have 
systematically larger FWHM(NB816) than ACS-detected LAEs, they may have a large UV 
continuum size and thus their surface brightness in Jgi4 falls below our detection limit. 

(6) UV continuum sizes of LAEs are not directly related to star formation properties 
such as Lja luminosity and Lya equivalent widths. 

We would like to thank both the Subaru and HST staff for their invaluable help, and 
all members of the COSMOS team. We would also like to thank the anonymous referee for 
his/her useful comments. This work was in part financially supported by JSPS (15340059 
and 17253001). 
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Table 1. SExtractor parameters for source detection 



Parameter 


Value 


Comment 


DETECT.TYPE 


CCD 


Detector type 


DETECT.MINAREA 


9 


Minimum number of pixels above threshold 


DETECT.THRESH 


1.6 


Detection Threshold in sigma 


ANALYSIS.THRESH 


2.0 


Limit for isophotal analysis in sigma 


FILTER 


Y 


Apply filter for detection 


T7TT 'TT~ T~> TVT A A JT? 

b 1L1 hjK-lN AMbj 


gauss_2 . 0_5x5 .conv 


Name of the filter for detection 


DEBLEND.NTHRESH 


32 


Number of deblending sub-thresholds 


DEBLEND.MINCONT 


0.015 


Minimum contrast parameter for deblending 


CLEAN 


Y 


Clean spurious detection 


CLEAN.PARAM 


1.0 


Cleaning efficiency 


MASK.TYPE 


CORRECT 


Correct flux for blended source 


PHOT.AUTOPARAMS 


2.5, 0.5 


MAG-AUTO parameters: Kron factor and minimum radius 


SATUR.LEVEL 


50000 


level (in ADUs) at which arises saturation 


MAG.ZEROPOINT 


25.936 


Magnitude zero-point 


MAG.GAMMA 


4.0 


gamma of emulsion (for photographic scans) 


GAIN 


1.0 


Detector gain in e _ /ADU 


PIXEL.SCALE 





Size of pixel in arcscc (0=use FITS WCS info) 


SEEING.FWHM 


0.11 


Stellar FWHM in arcscc 


BACK.SIZE 


64 


Background mesh size 


BACK.FILTERSIZE 


3 


Background filter size 


BACKPHOTO.TYPE 


GLOBAL 


Photometry background subtraction type 



Table 2. COSMOS z « 5.7 LAE sample 



LAE sample Number of LAEs Spectroscopic confirmation 



Total (Murayama 2007) 


119 


24 


Out of the ACS/F814W field 


34 


4 


In the ACS/F814W field 


85 


20 


ACS/F814W not-detected 


38 


3 


ACS/F814W detected 


47 


17 


Double 


2 


2 
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Table 3. ACS F814W Properties for the LAEs at z « 5.7 



Number 3, 


T b 


-Rhl 


Ism (estimated) d /?hl 


(estimated) d 




(mag) 


(arcsec) 


(mag) 


(arcsec) 


Detected 


2 


26.9 ± 0.1 


0.11 


26.818;$ 


10+ - 10 

u - iu -0.05 


4 


26.5 ± 0.1 


0.20 


26.ll°3 


1 n Qn 

0.35183° 


5 


26.7 ± 0.1 


0.13 


26.61°,^ 


o.i3i83° 


7 


27.5 ± 0.2 


0.10 


97 0+0.4 

z ' - o -0.6 


0-08l°3« 


13 


26.8 ± 0.1 


0.09 


26.7l°| 


o.o8i°:8« 


14 


25.7 ± 0.1 


0.24 


25.418;! 


0.381°;^ 


15 


26.0 ± 0.1 


0.24 


25.9181 


0.25l°3° 


20 


26.4 ± 0.1 


0.23 


26.1181 


0.381°:^ 


23 


27.1 ± 0.2 


0.13 


27.0133 


o.i3i°:^ 


30 


27.0 ± 0.1 


0.07 


26.9l8i 


0-081^1 


34 


26.4 ± 0.1 


0.25 


26.218J 


0.28l°3° 


35 


26.3 ± 0.10 


0.20 


26.2l8i 


n 1 s+o-13 

U - i8 -0.08 


39 


26.8 ± 0.1 


0.10 


26.718^ 


in+ 08 


40 


25.4 ± 0.1 


0.36 


q+0.6 


0.431™ 


41 


26.3 ± 0.1 


0.19 


26.3183 


o.i8l83^ 


42 


27.7 ± 0.1 


0.08 


27.5183 


n c i + °- 05 

u - ud -0.03 


43 


27.3 ± 0.2 


0.12 


27.018:^ 


n oq+0.40 
u -°°-0.20 


44 


27.1 ± 0.2 


0.14 


26.9181 


o.i8i°38 


45 


26.6 ± 0.1 


0.18 


26.418J 


0.28l°3| 


47 


26.1 ± 0.1 


0.16 


26.0+81 


U - iO -0.05 


49 


26.5 ± 0.2 


0.20 


26.0+83 


0.38+81 


50 


26.4 ± 0.1 


0.14 


26.4183 


n iq+0.08 

u - io -0.08 


51 


26.9 ± 0.2 


0.11 


26.9183 


o.i3l83^ 


55 


27.1 ± 0.2 


0.15 


26.5!^ 


0.25181 


60 c 


26.8 f 


0.94s 






60a 


26.9 ± 0.2 


0.15 


26.7183 


ri 9O+0.35 
u - zo -0.13 


60b 


27.3 ± 0.2 


0.11 


97 1 +0.5 

z ' • i -0.6 


0.10183°, 


68 


25.9 ± 0.1 


0.40 


25.7183 


0.35181 


69 


27.8 ± 0.2 


0.08 


27.5183 


0.031831 
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Table 3 — Continued 



Number a 


T b 
J 814 




i 814 (estimated) Iihl 


(estimated) 




(mag) 


(arcsec) 


(mag) 


(arcsec) 


71 


26 1 +01 


10 


26 l +0 - 2 

-L— 0.1 


08+ - 05 


73 


26 8 + 1 


1 3 


26 7 +0 - 4 


13+ 013 


74 


2fi q + n 2 


16 


26 6 +1 - 2 

ZjU .vj i 9 


35+ - 40 

U.Od_ 23 


75 


26 + 1 


1Q 


26 +0 - 3 


15+ - 10 


76 


26.7 ± 0.1 


0.10 


26.7+°' 3 

" —0.3 


0.08+S-Si 

" •"*-' — 0.05 


77 


27 7 + 1 

_ / . ( _1_ L/ . J. 


07 


27 5 +0 - 4 


05 +0 - 05 


78 


25.7 ± 0.1 


0.30 


25 4 + °i 


35+^ 

w.Otj_o 13 


81 


27 1 +01 


16 


26 9 +0 - 7 

zu - y -0.9 


28 +0 - 43 

u - zo -0.18 


83 




n DQ 

KJ.KJcr 


97 n+0.4 
Z ' - Z -0.4 


08 +0 - 13 

u - u °-0.05 


84 


26.2 ± 0.1 


0.15 


26.2 + no 
" — 0.3 


0.13+°^ 

—0.05 


85 


26 5 + 1 


1 7 


26 4 +0 - 3 

z,u -^-0.3 


15+ a10 

u - 1J -0.05 


8fi 


Zj 1 . VJ _l_ W . Zj 


n 1? 

W. -LZj 


97 n+0.4 
z ' - o -0.9 


o m+ - 28 

u - iu -0.05 


95 


27.0 ± 0.1 


0.11 


26.9 + ^ 4 

""-0.4 


o.io + °i° 

" "-0.05 


96 


26.6 ± 0.1 


0.19 


26.41°;^ 


2 c i +0 - 30 

U - ZJ -0.13 


98 


26.0 ± 0.1 


0.14 


26.0±g;l 


in+ - 08 

u - iu -0.05 


99 


26.8 ± 0.1 


0.16 


26.6^i 


o.is±g:Sg 


104 


26.5 ± 0.2 


0.20 


26.3 + °^ 


o-3o±g:?i 


105 


25.6 ± 0.1 


0.10 


25.618:1 


O.OSl™ 3 


107 


27.0 ± 0.2 


0.13 


26.818:2 


0.201^ 


110° 


26.4 f 


0.36 g 






110a 


26.5 ± 0.1 


0.17 


26.5l '; 4 


0.1818:13 


110b 


27.1 ± 0.2 


0.18 


26.8^; 9 7 


0.2818;*! 


Not-Detected 


1 


> 27.4 








6 


> 27.1 








9 


> 27.3 








16 


> 27.4 








19 


> 27.2 








21 


> 27.4 








24 


> 27.4 









Table 3 — Continued 



Number a / 8 i 4 b -Rhl c hu (estimated) d Rkl (estimated) d 
(mag) (arcsec) (mag) (arcsec) 



25 > 27.5 

26 > 27.4 

27 > 27.3 

28 > 27.3 

29 > 27.4 
38 > 27.3 
48 > 27.5 
54 > 27.5 

57 > 27.1 

58 > 27.4 

59 > 27.8 

61 > 27.5 

62 > 27.4 

79 > 27.6 

80 > 26.9 

87 > 27.2 

88 > 27.5 

89 > 27.3 

91 > 27.4 

92 > 27.6 

93 > 27.3 
97 > 27.2 

100 > 27.3 

101 > 27.6 

102 > 27.3 

103 > 26.9 
106 > 26.7 
108 > 27.4 

111 > 27.2 

112 > 27.3 
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Table 3 — Continued 

Number a 7 8 i 4 b -Rhl c hu (estimated) d Rul (estimated) d 
(mag) (arcsec) (mag) (arcsec) 

114 > 27.1 

a The LAE IDs given in Murayama et al. (2007). 

b AB magnitude. Lower-limits represent 3a significance with assum- 
ing a 1-arcsec diameter aperture. 

c Effective radius. 

d Estimated values by Monte Carlo calculations. 
e Double ACS sources. 

f Magnitude corresponding to the sum of flux of each source. 
Reparation between the double sources. 



Table 4. Photometric properties from Subaru Suprime-Cam Images for the LAEs at z « 5.7 



Number a 


i' h 


z' h 


iV5816 b 


FWHM(NB816) C 


e(NB816) d 


P^(NB816) 


L(Lya) 


EW (Lya) c 


z f 




(03"O) 


(0370) 


(03"O) 


(arcsec) 




(degree) 


(10 42 ergs s" 1 ) 


(A) 




Detected 


o 

z 


9fi ^ 
ZU.O 


9^ Q 


9/1 


i i 
i . i 


n 1 7 

U. 1 / 


Q 


s n 

o.U 


^JrU 




4: 


9zL ■=» 


9zL R 


93 9 
zo.z 


1 4 
i 


n ns 

u.uo 


1 R8 

1UO 


^i n 

Ol.U 


U i 




■J 


97 1 

Z 1 . 1 


9fi 9 
zu.z 


9zL zL 


1 9 

l.Z 


n D4 

u.u^ 


1 41 
i^i 


in 7 

1U. 1 


ou 




7 


9fi 


9^ 
zo.u 


9zL 


1 3 


nil 


1 7Q 
i / y 


7 Q 


3zL 
o^± 




lo 


97 
Z / .U 


QQ n 


93 Q 

zo.y 


1 1 
1 . 1 


n n7 


Do 


1 8 zL 

10.4: 


>> 1 7^ 

^> I/O 




1 4 
i*± 


9zL Q 
z^i. y 


9^ n 
zo.u 


93 
zo.u 


1 4 


U. lO 


74 


91 9 


oo 




1 ^ 
10 


9fi 9 
zu. z 


9^ 1 

ZO. 1 


9zL 


1 n 

l.U 


n nfi 


Q9 


7 


1 Q 
iy 




9D 
zu 


9fi ^ 
zu. o 


9fi 
zu.u 


94 7 


1 n 

i.u 


u.uo 


1UU 


s n 

O.U 


O 1 




93 


9^ zL 


9^ Q 
zu.y 


9zL zL 


1 4 
i .^t 


n ^7 

u.o / 


Q1 

c/l 


i n i 

1U. 1 


oo 




^n 


9^ Q 
zo. y 


9^ S 
zo.o 


9zL 3 
z^.o 


1 n 

l.U 


nil 


1 97 


19 9 


uu 




1 A 
34 


26.4 


or r* 

25.5 


24.7 


1.7 


0.31 


o 

8 


7.3 


on 

29 


5.681 


35 


25.5 


25.6 


24.4 


1.5 


0.41 


14 


9.6 


42 


5.663 


39 


25.9 


26.0 


24.2 


1.1 


0.16 


154 


13.6 


89 


5.718 


40 


25.5 


24.4 


23.9 


1.5 


0.09 


174 


13.3 


20 


5.690 


41 


27.1 


99.0 


24.5 


1.0 


0.16 


6 


11.2 


> 103 




42 


27.0 


27.3 


24.8 


1.1 


0.06 


137 


8.0 


> 72 


5.681 


43 


26.3 


26.1 


24.2 


1.3 


0.10 


166 


13.2 


91 


5.657 


44 


26.3 


26.0 


24.3 


1.4 


0.31 


67 


11.6 


77 


5.711 


45 


25.7 


25.9 


24.5 


1.4 


0.20 


97 


9.6 


55 


5.668 


47 


25.5 


25.2 


24.1 


1.4 


0.18 


128 


13.5 


42 


5.714 


49 


26.8 


28.3 


24.8 


1.1 


0.24 


176 


8.0 


> 84 




50 


26.0 


99.0 


23.5 


1.2 


0.06 


160 


27.9 


> 260 




51 


27.0 


99.0 


24.9 


1.2 


0.26 


22 


9.5 


> 88 





Table 4 — Continued 



Number a 


(<£3"0) 


z /b 
(0370) 


iVS816 b 
(03"O) 


FWHM(NB816) C 
(arcsec) 


e(NB816) d 


PA(NB816) 
(degree) 


L(Lya) 
(10 42 ergs s" 1 ) 


EW (Lya) e 

/ A \ 

(A) 


f 


00 


97 3 


9fi fi 
ZO.O 


94 ^ 


1 3 

1.0 


n in 
u. iy 


1 30. 

ioy 


Q 8 

y.o 


^ Q7 




fin 


9^ 4 


9^ ^ 

zo.o 


94 1 

Z^t. 1 


1 ^ 

l.O 


n 99 
u.zz 


7Q 
i y 


13 7 
lO. / 


^fi 

OU 


^ fifil 

O.UU1 


fi8 


9fi 4 


9fi 7 


94 ^ 


1 n 

i.u 


n nfi 

u.uu 


t o 


i n 4 


^111 

s> 111 


^ fifin 

o.uuu 


fiQ 


97 8 
z / .0 


qq n 

yy.u 


94 fi 


1 3 

l.O 


n in 


4Q 
iy 


1 n n 

1U.U 


>> i ni 

^> 1U1 




71 


9^ 4 


9^ 3 
zo.o 
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Table 4 — Continued 



Number a 


i' h 




NB8l6 b 


FWHM(NB816) C 


e(NB816) d 


PA(NB816) 


L(Lya) W (Lya) c z sp f 
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118 


25.6 


25.5 
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0.17 


53 


22.8 94 


119 


25.8 


25.2 


24.5 
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0.36 


41 


8.2 25 



a The LAE IDs given in Murayama et al. (2007). 

b AB magnitude. An entry of "99.0" indicates that no excess flux was measured. 
c Spatial sizes measured on NB816 images. 
d Ellipticity measured on NB816 images. 
e Lower-limits represent la significance. 
Spectroscopic redshift. 
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Table 5. Symbols and references for Figure 18 



Symbol 


Description 




References 


Large red filled circle 


COSMOS LAEs atz« 5.7 




this study 


Small black filled circle 


i band dropout galaxies in UDF, 
GOODS-N, and GOODS-S 


UDF-P, 


Bouwens et al. (2006) 


Small red filled circle 


Spectroscopic confirmed LAEs among 
dropout galaxies in GOODS-S 


i band 


Bouwens et al. (2006) 
ESO GOODS-S database 11 


Small blue filled circle 


Spectroscopic confirmed LAEs among 
dropout galaxies in GOODS-S 


% band 


Bouwens et al. (2006) 
ESO GOODS-S database 11 


Red plus 


LAE at z = 5.7 




Bunker et al. (2003) 


Black asterisk 


% band dropout galaxies in UDF 




Bunker et al. (2004) 


Red open circle 


LAEs at z ~ 5.8 




Stanway et al. (2004a) 


Red cross 


LAEs in GOODS-N 




Stanway et al. (2004b) 


Blue open square 


Non emission-line galaxies at z ~ 6 




Dow-Hygelund et al. (2007) 


Red open square 


Spectroscopic confirmed LAEs at z ~ 6 




Dow-Hygelund et al. (2007) 



Spectroscopic properties for the sample from Bouwens et al. (2006) were taken from the ESO 
GOODS/VIRMOS spectroscopic database at http://www.eso.org/science/goods/ 
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150.5 150 149.5 

RA (deg.) 

Fig. 1. — Spatial distribution of our sample of 119 LAEs (M07). The whole area is 1.95 square 
degree that is mapped with our Suprime-Cam observations. The ACS image is available for 
1.64 square degree, its footprint is indicated by the solid gray line. Masked out regions are 
shown by filled gray circles or thin gray lines. The 47 LAEs detected with ACS are shown 
by large filled circles while 38 LAEs undetected in the ACS images are shown by crosses. 
The remaining 34 LAEs shown by dots fall outside the HST/ACS field. 
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Fig. 3. — Thumbnails of all 85 LAEs analyzed in this paper. Each panel has a size of 5"x 
5". Red ellipses overplotted on the smoothed ACS, NB816, and z' images are half light 
ellipses of the detected LAE components in the ACS image. Green and blue ellipses are 
ACS sources excluded from the sample by eye inspection and by rejection of foreground 
neighbors, respectively. 
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Fig. 4. — Distribution of the half-light radius in the ACS images, Rrl, for the 47 LAEs 
detected with ACS by a dotted histogram. LAEs with spectroscopic confirmation are shown 
by a solid histogram. 

The PSF size (0.11 arcsec) derived from stars is shown by the dashed line. 
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Fig. 5. — Frequency distribution of the ACS magnitude 7 814 by a dotted histogram. LAEs 
confirmed by our follow-up spectroscopy are shown as a solid histogram. 
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Fig. 6. — Frequency distribution of the ACS magnitude I$i4 for the not-detected LAEs 
(dotted histogram); 3a upper limits are given for each LAE. LAEs confirmed by our follow- 
up spectroscopy are shown by a solid histogram. 
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7 814 (magauto) 

Fig. 7. — Diagram between half-light radius -Rhl and the ACS magnitude Isu- The LAEs 
confirmed by our follow-up spectroscopy are shown by filled circle. For double-component 
LAEs, each component is plotted. The 50% detection completeness for exponential disk 
objects estimated by Monte Carlo simulation is indicated by the dashed curve. 



-61 - 



i — i — i — r 



i — i — i — r 

m 



i — i — i — r 

Detected 



10 — 



5 — 




10 



I I I I 



r 



h" 1 p i i i i i 







Not detected _ 



ft 



i" 
i 



r-l 
i i 



J I L 



iLtL 



I I ■ 



m □ 



J I l_L_L 







FWHM(NB816) (arcsec) 



Fig. 8. — Frequency distributions of size in the narrow-band filter NB816 FWHM(NB816) 
for the detected and not-detected LAEs by dotted histograms. LAEs confirmed by our 
follow-up spectroscopy are shown by solid histograms. 
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Fig. 9. — Diagram between size in the narrow-band filter NB816, FWHM(NB816) and the 
half-light radius in the ACS image, /?hl- LAEs confirmed by our follow-up spectroscopy 
are shown by filled circles. For double-component LAEs Rul is plotted for each component 
while the same FWH M (NB816) is adopted for the two components as they are not resolved 
in the NB816 images. 
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Fig. 10. — Results of the Monte-Carlo simulation: the relation between estimated total 
magnitudes by the simulation and measured magnitudes of our LAE sample (upper panel) 
and the relation between estimated half light radii by the simulation and measured half light 
radii (upper panel). 
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Fig. 11. — Transmission curves for the filters used in our analysis; F814W (blue) for 
ACS/HST and i' (green) , NB816 (magenta), and z' (red) for Suprime-Cam/Subaru. The 
CCD sensitivity is taken into account for each filter transmission curve. Model spectrum 
(black) of a LAE at z ~ 5.7 with a rest-frame EW of 30A is also plotted. 
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NBQ16 (3"0) 



Fig. 12. — Diagram between 7 814 and NB81Q. LAEs confirmed by our follow-up spectroscopy 
are shown by filled circles. For the double-component LAEs, Isu is the total magnitude of 
the two components. Their A^B816 magnitudes are measured with a 3-arcsec aperture that 
included the two components. 
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Fig. 13. — Diagrams between J 814 and %' (left), z' (right). LAEs confirmed by our follow-up 
spectroscopy are shown by filled circles. For the double-component LAEs, J 8 i4 is the total 
magnitude of the two components. Their i' and z' magnitudes are measured with a 3-arcsec 
aperture that included the two components. 
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Fig. 14. — Diagrams between the corrected z(I S u) and z' . LAEs confirmed by our follow-up 
spectroscopy are shown by filled circles. For z(I 8U ) of the double-component LAEs, their 
total magnitudes of the two components are used. Their Subaru z' magnitudes are measured 
with a 3-arcsec aperture that included bot components. 
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Fig. 15. — Diagrams between /?hl and L(Lya). LAEs confirmed by our follow-up spec- 
troscopy are shown by filled circles. For double-component LAEs Rul is plotted for each 
component while the same L(Lya) is adopted for the two components, because they are not 
resolved in the ground-based images. 
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Fig. 16. — Diagrams between i? H L and EW (Lya). LAEs confirmed by our follow-up spec- 
troscopy are shown by filled circles. LAEs confirmed by our follow-up spectroscopy are shown 
by filled circles. For double-component LAEs -Rhl is plotted for each component while the 
same EW(Lja) is adopted for the two components, because they are not resolved in the 
ground-based images. 
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Fig. 17. — Composite (stacked) ACS image of all LAEs. The surface brightness profiles 
(thick lines) , the model profiles (red lines) and residuals (thin likes) are plotted. The light 
profile of the PSF (dashed line) derived from stars in the ACS image is also indicated for 
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Fig. 18. — Diagrams between /?hl and Z850 for z ~ 6. The left panel is the same as the 
right panel except that our data is excluded for clarity. Spectroscopically confirmed LAEs 
are shown by the red symbols and spectroscopically confirmed LBGs without Lya emission 
are shown by the blue symbols. Our LAEs are shown by large red filled circles. For the 
double-component LAEs, each component is plotted. See Table 5 for explanation of the 
remaining symbols. 



